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FIELD OF THE INVENTION 



The present invention relates to a method for measuring the wettability of rock 
samples by measuring the relaxation time by low-field nuclear magnetic resonance. 

The method finds applications notably for analysis of rocks taken fi-om an 
imderground formation containing or likely to contain hydrocarbons. 

Knowledge of various parameters, and notably of the wettability of the rocks, is 
also usefiil notably for enhanced recovery of a formation, by injection of a fluid imder 
pressure, and when the fluid (liquid or gas) best suited for effluents displacement is to 
be determined by means of preliminary tests. 

The invention also finds applications in civil engineering for formation hydrology 
in order to evaluate the degree of pollution of formations for exanq>le, or in the building 
trade for testing building materials notably in order to determine waterproofing 
treatments for example. 



There are well-known methods for determining the wettability of rocks to the water 
and oil that they may contain, comprising carrying out rock drainage cycles, ie. 
displacement of fluids intended to decrease the water saturation, followed by imbibition, 
this term relating to a displacement of the fluids allov^g to increase the water 
saturation (Sw) of the rock. The capillary pressure Pc at one point is defined as the 
difference at equilibriiun between the pressure P(oil) of the oil and the pressure P(water) 
of the water. This parameter makes sense only if the two fluids are in the continuous 
phase in the porous medimn. For a water wet medium, only the positive values make 
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sense. On the other hand, vAien the medium has a mixed wettability, the fluids can 
remain in the continuous phase for the positive as well as for the negative capillary 
pressures (Pc). 

For an application of this type, a complete capillary pressure measuring cycle 
5 generally comprises (Fig.l) : 

a) positive primary drainage of an initially 100 % water-saturated sample (curve 1) ; 

b) positive imbibition (curve 2) ; 

c) negative imbibition (curve 3) ; 

d) negative drainage (curve 4) ; and 

10 e) positive secondary drainage (curve 5). 

There are various types of devices allowing the curves of Fig.l to be drawn. 

In a device referred to as « porous plate » device, notably described in patent US- 
4,506,542, the porous rock sample containing two fluids in the continuous phase is 
placed in an elongate cell ended at its two opposite ends by capillary barriers permeable 
15 to a first fluid. This first fluid is injected under pressure through the first membrane and 
the pressure difference between the injection pressure and the pressure of the fluid 
discharged at the other end is measured. The pressures of the two fluids and the 
capillary pressure Pc are constant all along the sample, and the saturation is assumed to 
be uniform. 

20 It is also known to carry out progressive-speed centrifugation by means of 

centrifugation devices such as those described for example in patents FR-2,772,477 
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(US-6,1 85,985) or FR-2,763,690, or in patents EP-603,040 (US-5,463,894) or FR- 
2,798,734 filed by the applicant. 

The sample initially saturated with a first Uquid (brine for example) is placed 
(Fig. 10) in an elongate vessel or cup A containing a second fluid of different density 
5 (oil for example). Cup A is fastened to the end of a rotating arm B and a centrifiigal 
force is applied thereto so as to study the displacements of the fluids in the sample 
during at least two distinct stages. During a first drainage stage, the assembly is then 
subjected to a centrifiigal force applied along the length of the vessel so as to exert on it 
an expulsion force which tends to cause part of the first fluid to flow out. The pressure 

10 field created by centrifiigation is expressed as a fimction of the density p, of the radius 
R and of the angular velocity w, by the relation : V^ w^** . p (Rmax^ - R^), for each fluid. 
We impose that the pressure of the two fluids at the sample outlet is the same and is 
equal to zero at the outlet. At the same time, the second fluid flows into the sample. The 
two fluids move in the sample until a position of equihbrium is reached where the force 

15 due to the capillary pressure in the pores compensates the centrifiigal force exerted. A 
measuring sonde is placed in the cup on the side of the sample. The sonde can be of 
capacitive type for example and coiiq>rise a metal rod insulated by a ceramic sheath. 
The capacity between the rod and the conducting fluid (brine), which is proportional to 
the immersed height, is measured. With this measuring means, the measuring accuracy 

20 is 1.5 % of the pore volume. The sonde detects the position of the interface between the 
two liquids in the cup and transmits the measuring signals to a measuring signal 
acquisition and control automaton E includiag hydraulic liquid circulation means and an 
acquisition device. 
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During the re-imbibition stage, the velocity is decreased so as to study the return of 
the initial fluid into the sample. The local saturations measured with this type of device 
are calculated by an inversion programme from the total amount of water expelled from 
the sample. 

5 According to another method, referred to as « dynamic » method, a sample is placed 

in an elongate cell ended at its two ends by water-permeable membranes. At a &st end, 
oil under pressure is directly injected into the enclosure. Water is also injected, but this 
injection is carried out through the membrane and at a lower pressure. At the opposite 
end, the oil is directly discharged \^ereas the water flows out through the terminal 

10 membrane. Adjustment of the oil and water injection rates allows the capillary pressure 
to be the same at the inlet and at the outlet of the enclosure, which leads to a uniform 
saturation that can be deduced from the fluids balance. The capillary pressure is 
obtained for example by measuring the difference between the pressure of the oil and of 
the water at the enclosure outlet. Such a method is notably described by Brown H. W. in 

15 « Capillary Pressure Investigations », Petroleum Transaction AIME, vol. 192, 1951. 
Examples of implementation are for example described in patents EP-729,022 (US- 
5,698,772) or EP-A-974,839 (US-6,229,312) filed by the appUcant. 

A method referred to as semi-dynamic method is also known, wherein a rock 
sample imbibed with a first fluid is confined in a closed cell, another fluid under 
20 pressure is injected at a first end of the enclosure and the opposite end is swept by a 
low-pressure fluid circulated by pmnpiag means which carries the drained fluid outside. 
The device comprises means for measuring the pressure and the saturation of the 
sample, the amount of fluid discharged and the electric resistivity of the sample. This 
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method is implemented for example in patent FR-2,708,742 (US-5,679,885) filed by the 
applicant. 

These drainage and imbibition curves being established, it is well-known to 
calculate the wettability indices WI from the surface areas Ad and Ai delimited by the 
5 positive and negative capillary pressure curves, as shown in Fig.l, by the relation 

SUMMARY OF THE INVENTION 

The method for measuring the wettability of a porous rock sample in the presence 
of water and oil according to the invention comprises determining the water wet pore 
10 surfece and the oil wet pore surfece when the sample is saturated with water and oil, and 
calculating the wettability index by combination of the values obtained for said 
sur&ces. 

Determination of the water wet pore sur&ce and of the oil wet pore sur&ce when 
the sample is satxarated with water and oil is obtained for example by means of 
15 measurements of relaxation times of the sample placed in a nuclear magnetic resonance 
device. 

The wettability index is for example determined by the relation : 



= ca/ ^ cii/ relation = log.o — — 



where SM^ is the water wet pore surfece and SMo is the oil wet pore surfece when the 
porous medium is saturated with water and oil. 
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According to an implementation mode, the wettability index is deteraiined by the 
following operations : 

a) measuring the characteristic relaxation times of the water-saturated sample ; 

b) measuring the characteristic relaxation times of the oil in the sample in the presence 
of water, in a zone close to saturation (Swir) of the sample ; 

c) measuring the characteristic relaxation times of the water in the sample in the 
presence of oil, in a zone close to residual saturation (Sor) ; 

d) measuring the relaxation times of the sample in a state where its 100 % oil 
saturation point is reached ; and 

e) combining the measurements of the various relaxation times so as to obtain said 
index. 

According to a preferred implementation mode, the relaxation times of stages a) to 
c) are determined after subjectiag the sample to centrifiigation. 

According to a preferred implementation mode, the relaxation times of stage d) are 
determined after forced displacement of the fluids in the sample placed in a containment 
ceU. 

An oil whose intrinsic relaxation time (Tb) is as great as possible and as close as 
possible to that of water, dodecane for example, is preferably selected. 

The characteristic relaxation times are for exanq>le those corrresponding to either 
the saturation curves maxima, or to mean values of said curves. 
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Experience shows that the measurement of Inmr obtained by means of the method is 
as sensitive but that it requires much less time and that it is applicable to a large nxmiber 
of samples. 

BRIEF DESCRIPTION OF THE HGURES 

5 Other features and advantages of the method and of the device according to the 

invention will be clear from reading the description hereafter of an embodiment given 
by way of non limitative example, with reference to the accompanying drawings 
wherein : 

Figure 1 shows examples of capillary pressure curves necessary for determination of 
10 the conventional wettability index Iusbm (United States Bureau of Miaes), 

- Figures 2a and 2b respectively show the distribution of relaxation times Ta for a 
water wet rock saturated with water and oil, and a representation of the phase 
distribution (matrix in hatched lines, water in light grey and oil in darker grey), 

- Figures 3a and 3b respectively show the distribution of relaxation times T2 for a rock 
15 of intermediate wettabiUty saturated with water and oil (same central part as in Figure 

2), and a representation of the phase distribution (matrix in hatched lines, water in light 
grey and oil in darker grey), 

- Figure 4 shows, on capillary pressure curves, the saturation values where the NMR 
relaxation measurements are performed, 

20 - Figiue 5 shows an example of distribution of relaxation times T2 at the saturation 
values given in Figure 4 (from the bottom up : 100 % water, Swi, Sor and 100 % 
dodecane). 
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- Figure 6 shows the relation between index and the standard index Iusbm^ 

- Figure 7 diagranunatically shows a containment cell with a permanent heating 
circuit by fluid circulation, and 

- Figure 8 diagrammatically shows the layout of a centrifiigation device used for 
5 implementing the method. 

DETAILED DESCRIPTION 

It may be reminded that the NMR analysis technique essentially consists in 
applying to an object to be tested a first static magnetic polarization field Bo intended to 
align the initially randomly oriented hydrogen protons nuclei in the direction of the 

10 field, then a second impulsive magnetic field oscillating at the Larmor firequency, 
perpendicular to the first one, created by coils excited by a control signal to carry out a 
nuclear magnetic resonance experiment. When this inapulsive field stops, the retum of 
the nuclei to their initial state or relaxation generates electromagnetic signals (echoes) 
which are detected and analyzed. The presence of physical parameters of the object is 

15 determined firom the amplitude characteristics of these signals. 

We propose defining a new wettability index constructed by combming values of 
the water wet pore surface SMw and of the oil wet pore surfece SMo when the porous 
medium is saturated with water and oiL The index can be calculated for exanq>le by the 
relation as follows : 

j^=^M.z^ (1.) 

SM„ + SM„ 
20 or by the relation : 
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Quantities SM^ and SMo are obtained for example by measuring the dominant 
relaxation time in a low-field nuclear magnetic resonance experiment. 

In fact, in such an eiq^eriment, the fundamental relation connecting the longitudinal 
Ti or transverse T2 relaxation time to surfece S and to the voliraie of this pore is as 
5 follows : 



^ 1.2 ^ \B,2B ^ 

where Tib^b is the relaxation time of the fluid outside the porous medium. Basically, this 
relation comes from the feet that the molecules in the pore imdergo difiiision motions 
and interact with the porous wall during the measuring time (the typical relaxation time 
is of the order of 100 ms). The surfece interactions are designated by coefficient pi^ 

10 referred to as surfece relaxivity. Besides, we have disregarded a gradient term that is 
important when measurements are performed at a high magnetic field (> 0. IT). Relation 
2 is strictly valid for a pore saturated with a single fluid. In general, natural porous 
media have a pore size and therefore ratio SA^ distribution. A relaxation times 
distribution is thus generally observed, but this does not modify the method described 

15 here. 

When two fluids are present in a pore within the porous medium, the same physical 
difiiision mechanism is valid, but the distribution of the two fluids in relation to the pore 
surfece is of crucial importance. For example, when the medium is water wet, the water 
is at the surfece and interacts therewith, whereas the oil is inside the pore and does not 
20 interact with the surface. When the distribution of the relaxation times is measured for 
such a system, the result of Figure 2a is observed. The water has shorter relaxation 
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times than in the case of a 100 % water saturation, because the volume of water V has 
decreased (see equation 1). The oil is characterized by a relaxation time that is the same 
as when the measurement is performed outside the porous medium (Tib^) because there 
is no interaction with the pore surface. Besides, it is well-known that interactions at the 

5 water/oil interface do not lead to a relaxation that is different from that of the oil outside 
the porous medium, provided that the oil considered is light (refined oil such as 
dodecane). Diagrammatically, the phase distribution is shown in Figure 2b. When the 
same porous medixmi has a different wettability and is not clearly water wet, the 
distribution of the relaxation times is modified (Figure 3a). Essentially, the 

10 characteristic time of the oil will be shorter as a result of an interaction of the oil at the 
pore sur&ce because the water is no longer the closer to the surface. The possible 
distribution of the two fluids is shown in Figure 3b. Such an effect has been obtained 
using standard procedures allowing to reproduce the surfece properties of rocks from a 
petroleum reservoir : cleaning, then aging in the reservoir oil at the temperature of the 

15 reservoir at the irreducible water saturation. 

We thus generalize relation (2) to a two-phase water-oil system. For the water, we 
have : 



1 1 SM^ 



and for the oil : 



1 1 SM„ ... 



We can thus determine quantities SMw and SMo by measuring relaxation times Ti or 
20 T2 in the porous medium, relaxation times T,b or T2B of the fluids outside the porous 



11 



mediiUD, liquid volumes Vo or Vw present in the porous medium. The interaction 
constants p can be determined separately, but we will see that only the ratio of these 
constants is necessary. 

If we introduce relations 3 and 4 in base relation 1, we obtain : 



Sw(^---) + CpSo(--—) 

5 where Sw=VwA^p is the water saturation, 
So=VoA^p is the oil saturation^ and 
Vp the pore volume. 

Cp = — is the ratio of the relaxivities for the water wet (w) and oil wet (o) 

Po 

surfaces. To simplify the notation, indices 1 and 2 have been suppressed. 

10 The most suitable saturations for measuring quantities SM^ and SMo thus remain to 

be determined. For reasons linked with the calculation of the relaxation tune 
distributions, we select the irreducible water saturation (Swi, Figure 4) to determine the 
oil wet surfece in the presence of water (SMw), and the residual oil saturation (Sor, 
Figure 4) to determine the water wet sur&ce in the presence of oil. We also select an oil 

15 whose intriosic relaxation time (Tb) is as great as possible and as close as possible to 
water. In feet, if Tb is too small, the method will be limited to porous media whose ratio 
SA^ is high (small pore sizes), or to media whose surface relaxivity is high. Dodecane 
for example is a refined oil which is suited to the proposed measurement because its 
intrinsic relaxation time (Tbo) is 1 s, close to the relaxation time of water (Tbw, about 2.7 

20 s). Crude petroleum oils should generally not be used because their relaxation times are 



12 



too small and they also have a relaxation time distribution which considerably hinders 
analysis. The practical details of the calculation are given below. 

From the distribution of the relaxation times at the different saturations (see the 
example of Figure 5), the relaxation time corresponding to the dominant peak (Tw at 
5 Sor, To at Swi) can be easily determined, and these values are used in relation 5. It is 
well-known that determination of the dominant relaxation time is robust and depends 
little on the distribution calculation process. It can be noted that the dominant peak 
corresponds to that of the oil at Swi and of the water at Sor. The presence of the second 
fluid at these different saturations is important from a physical point of view, but it 

10 disturbs calculation very little. We can show that whatever the relaxation time of the 
second fluid present in small amount (water at Swi, oil at Sor), the relaxation time of the 
dominant fluid does not fluctuate much. This aspect is importaut for the calculation 
robustness. Figure 5 shows the distribution of the transverse relaxation times T2. The 
longitudinal relaxation time Ti can also be used, but this measurement is less favourable 

15 in general because Ti > T2 but Tib = T2B for a given fluid. 

In order to detennine the surface relaxivity ratio Cp, we use the dominant relaxation 
times T^ioo and Tdoo respectively at the two saturations Sw=100 % and So=100 %. The 
formula used is as follows : 

" Po lfT„,^-VT^ 
Figure 5 illustrates the two measurements performed. Thus, in total, the relaxation 
20 times distributions have to be measured at four saturations : 100 % brine, 100 % oil and 
the extreme saturations Swi and Sor. We can therefore use the centrifugation technique 
to reach, starting from Sw=100%, Swi, then Sor. A single rotation stage at maximum 
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speed allows these values to be rapidly reached. For measurement at So=100%, we use 
a succession of miscible displacements by injecting an alcohol, then the refined oil- 

Implementation 

For implementing the method, we use for example an NMR measuring device 
5 conventionally con^jrising (Fig.7) magnets 1 with an air gap in which coils 2 are 
arranged. The coils are connected to an electronic box to generate the excitation signals 
(signals in the radiofrequency range for example) creating an oscillating electric field 
and to acquire the response of sample S to the excitation signals. A device of this type is 
used for example in patent appUcation FR-A-2,823,308. 

10 Practical implementation example 

The wettability index is obtained by carrying out for example the following 
succession of stages with a reservoir sample whose wettability is intermediate : 

1. The sample is 100% saturated with reservoir brine and placed in the NMR 
measuring device (as shown in Figs.7-8 for example) to measure the distributions of 

15 relaxation times Ti, T2 and to deduce parameter TwlOO of equation 6 ; 

2. It is then placed in the centrifugation device (as schematized in Fig. 10) and 
centrifiiged at maximum speed in a vessel filled with oil xmtil the irreducible water 
saturation Swir is reached ; the distributions of relaxation times Ti, T2 are then measured 
in the NMR device so as to deduce To of equation 5 ; 

20 3. After being placed again in the centrifugation device, the sample is centrifiiged 

at maximum speed imtil the residual oil saturation Sor is reached ; it is then transferred 
into the NMR device again to measure relaxation times Ti, T2 and to deduce T^ ; 



4. The sample is thereafter placed in a containment cell such as those described in 
the aforementioned patents and a water and oil miscible solvent, then oil (dodecane for 
example) is injected until the 100% oil saturation is reached. This point bemg reached, 
relaxation times Tj, T2 are again measured in the NMR device and the value of Toioo 
5 required for detennination of Cp is deduced. 

We thus know all the parameters for calculation of wettability coefficient Inmr. 
Comparison with index Iusbm 

Comparison between wettability index Wm and the new index Inmr shows a good 
correlation between these two quantities (Figure 6) for reservoir rocks. The 
10 measurement of Iusbm can therefore be advantageously replaced by a measurement of 
Inmr, which is as sensitive but requires much less time and is appUcable to a large 
number of samples. 

We have described an implementation example where a NMR type relaxometry 
method is used to measure the surfece and the volimtie of the pores. This method is 
15 however not limitative. Any other analysis method can be used, notably 
cryomicroscopy. 



